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ABSTRACT 

We select a sample of young passive galaxies from the Sloan Digital Sky Survey Data 
Release 7 in order to study the processes that quench star formation in the local universe. 
Quenched galaxies are identified based on the contribution of A-type stars to their observed 
(central) spectra and relative lack of ongoing star formation; we find that such systems ac- 
count for roughly 2.5 per cent of all galaxies with log (M/M ) > 9.5, and have a space density 
of ~2.2 X 10" 4 Mpc~ 3 . We show that quenched galaxies span a range of morphologies, but that 
visual classifications suggest they are predominantly early-type systems. Their visual early- 
type classification is supported by quantitative structural measurements (Sersic indices) that 
show a notable lack of disk-dominated galaxies, suggesting that any morphological transfor- 
mation associated with galaxies' transition from star-forming to passive — e.g. the formation 
of a stellar bulge — occurs contemporaneously with the decline of their star-formation activ- 
ity. We show that there is no clear excess of optical AGN in quenched galaxies, suggesting 
that: i) AGN feedback is not associated with the majority of quenched systems or ii) that the 
observability of quenched galaxies is such that the quenching phase in general outlives any 
associated nuclear activity. Comparison with classical post-starburst galaxies shows that both 
populations show similar signatures of bulge growth, and we suggest that the defining charac- 
teristic of post-starburst galaxies is the efficiency of their bulge growth rather than a particular 
formation mechanism. 
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1 INTRODUCTION 

Large surveys such as the Sloan Digital Sky Survey (SDSS; York 
et al. 2000) and 2dF Galaxy Redshift Survey (Colless et al. 2001) 
have placed tight constraints on properties of the local galaxy pop- 
ulation and, increasingly, spectroscopic surveys out to z ~ 1 (e.g. 
DEEP2 [Davis et al. 2003] and zCOSMOS [Lilly et al. 2007]) in- 
form our phenomenological picture of how galaxies have evolved 
over the past 8 billion years. For example, we know that the lo- 
cal galaxy population is bimodal in colour: passive galaxies form 
a tight 'red sequence' in the colour-stellar mass plane while star- 
forming galaxies populate a more dispersed 'blue cloud' (e.g. Strat- 
eva et al. 2001 ; Baldry et al. 2004), and these sequences are in place 
already by at least z = 1 (e.g. Cirasuolo et al. 2007; Cooper et al. 
2008; Williams et al. 2009; Brammer et al. 2009; Muzzin et al. 
2012). The number density of red-sequence galaxies is observed 
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to evolve significantly from z ~ 1 to the present (Bell et al. 2004; 
Faber et al. 2007), suggesting that star-formation quenching plays 
an important role in galaxy evolution over this epoch. 

However, it does not appear that simple evolutionary scenar- 
ios, e.g. gas exhaustion, can explain the transition of galaxies from 
star-forming to quiescent in the absence of additional physical pro- 
cesses. There is a correlation between star formation and galaxy 
structure such that passive galaxies are typically spheroid domi- 
nated, while star-forming galaxies are generally disky (e.g. Bam- 
ford et al. 2009; Skibba et al. 2009, but see also Masters et al. 2010), 
and both star formation and morphology are connected to galaxies' 
environment (e.g. Lewis et al. 2002; Goto et al. 2003; Kauffmann 
et al. 2004; Balogh et al. 2004, and references therein). The buildup 
of passive galaxies since z ~ 1 appears to be mirrored by a change 
in galaxies' light profiles (e.g. McGee et al. 2008, 201 1; Bell et al. 
2012), and therefore the evolution of both star-formation rate and 
morphology must occur in such a way as to jointly preserve the 
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relationships between galaxy structure, star formation and environ- 
ment from high redshift to the present. 

Numerous physical mechanisms have been proposed to ex- 
plain galaxies' transition from star-forming disks to passive 
spheroids with time, both in terms of their internal properties and 
their surrounding environments. Mergers are thought to play an im- 
portant role in both assembling stellar mass and triggering efficient 
star formation in groups and the field, where low encounter veloc- 
ities are favourable to strong tidal interactions. Major mergers (be- 
tween equal-mass galaxies) can rapidly transform disk-dominated 
galaxies into spheroids (Toomre & Toomre 1972; Negroponte & 
White 1983; Mcintosh et al. 2008, but see also Hopkins et al. 2009), 
while minor mergers promote bulge growth by exciting instabil- 
ities in otherwise stable disks. There is evidence that tidal inter- 
actions can drive radial gas flows in disks, leading to centrally- 
concentrated starbursts (Mihos & Heraquist 1994; Ellison et al. 
2010; Scudder et al. 2012). Such inflows can also fuel accretion 
onto a central supermassive black hole; subsequent energetic feed- 
back may be sufficient to prevent additional gas cooling (e.g. Kauff- 
mann & Haehnelt 2000; Croton et al. 2006; Hopkins et al. 2008b; 
van de Voort et al. 2011), efficiently smothering star formation. 
Environmental processes tied either to galaxies' interaction with 
the intra-cluster medium (i.e. ram-pressure or viscous stripping; 
Gunn & Gott 1972; Nulsen 1982) or the properties of their host 
dark-matter haloes (so-called "halo quenching"; Keres et al. 2005; 
Dekel & Birnboim 2006) may also play an important part in reg- 
ulating star-formation activity. Increasingly, a combination of ex- 
ternal mechanisms and internal, secular processes are also invoked 
in order to reproduce the observed properties of the galaxy pop- 
ulation. Disk instabilities can lead to the formation of a bar (e.g. 
Sellwood & Wilkinson 1993, and references therein) which can, in 
turn, drive radial gas flows and the growth of a secular (pseudo) 
bulge. Observations in the local universe suggest that this secular 
bulge growth is a nearly ubiquitous feature in disk-dominated sys- 
tems (Kormendy et al. 2010). 

Confronting our theories of galaxy evolution with observa- 
tional data is a critical step in developing a refined, physical un- 
derstanding of the galaxy population as we see it today. Whole- 
sale comparison of the observed and theoretical galaxy population 
has been significantly aided by the development of so-called semi- 
analytic models (e.g. Kauffmann, White & Guiderdoni 1993; Cole 
et al. 1994; Somerville & Primack 1999), which 'paint' galaxies 
onto high resolution A'-body simulations using analytic prescrip- 
tions to capture the behaviour of relevant physical processes. Such 
models enjoy success in reproducing the bulk properties of galax- 
ies at z ~ 0, such as the luminosity and stellar mass functions (e.g. 
Bower et al. 2006; De Lucia et al. 2007) and the colour-magnitude 
relation (e.g. Weinmann et al. 2006b; Henriques et al. 2009). How- 
ever, despite these successes, several details of the galaxy popula- 
tion are poorly reproduced. Galaxies assemble too quickly, leading 
to a significant over-prediction of the stellar mass function at z > 1 
(e.g. Somerville et al. 2008; Guo et al. 201 1) and an offset between 
the ages of simulated and observed galaxies, particularly at low 
mass (e.g. Fontanot et al. 2009; Weinmann et al. 2012). Environ- 
mental processes are too efficient, resulting in an overproduction of 
passive low-mass galaxies in high-density environments by up to 
50 per cent (e.g. Weinmann et al. 2006b; Kimm et al. 2009). More 
generally, models seem to under-predict star-formation rates by a 
factor of ~4 at z < 1 (e.g. Daddi et al. 2007), suggesting our under- 
standing of mechanisms that regulate star formation is incomplete. 

One way to further constrain formation of the passive pop- 
ulation is to study directly the properties of galaxies in transition 



between star-forming and quiescent; however, such galaxies are 
expected to be rare in the local Universe. Nevertheless, samples 
of so-called "K+A", "E+A" or "post-starburst" galaxies' can be 
identified based on the presence of strong Balmer-line absorption — 
characteristic of A-type stars — and a lack of optical emission lines, 
and have recent star-formation histories characterised by intense 
star-formation followed by a precipitous decline. While such galax- 
ies were first studied in the context of clusters (e.g. Couch & 
Sharpies 1987), they have since been shown to inhabit a broad 
range of environments (Zabludoff et al. 1996; Blake et al. 2004; 
Hogg et al. 2006; Nolan, Raychaudhury & Kaban 2007), support- 
ing the idea that their formation is not necessarily governed by 
cluster-specific processes. Indeed, photometric observations have 
shown that anywhere from 30 to 80 percent of post-starburst galax- 
ies are morphologically disturbed (Yang et al. 2004b; Goto 2005; 
Pracy et al. 2009), highlighting mergers as one of the likely drivers 
of post-starburst galaxy formation. However, the high burst-mass 
fractions inferred for some post-starburst galaxies, anywhere from 
20-60 per cent (e.g. Kaviraj et al. 2007), suggest that they be biased 
towards relatively extreme formation scenarios, and may not reflect 
the evolution of "average" galaxies. 

In this paper we exploit a wealth of SDSS data to identify 
a large sample of galaxies in transition between star forming and 
passive, hereafter referred to as 'quenched', facilitating a direct 
study of their properties. This work builds upon the substantial 
body of literature characterising the integrated properties of the 
passive galaxy population (e.g. colour or star-formation rate; Wein- 
mann et al. 2006a; van den Bosch et al. 2008; Wetzel, Tinker & 
Conroy 2012) by focusing on a subsample of passive galaxies that 
have only recently stopped forming stars. In studying the frequency, 
physical properties and environments of such systems our goal is 
to constrain the processes that lead to their formation and, more 
generally, the physical mechanisms that contribute to the buildup 
of the red (passive) sequence over time. In this paper we iden- 
tify quenched galaxies based on presence of relatively young stel- 
lar populations coupled with low instantaneous star-formation rates 
(SFRs), and show that quenched systems bear a close resemblance 
to the population of galaxies already on the red sequence. In a forth- 
coming paper (Mendel et al. in prep.) we will discuss more directly 
the relationship between quenched galaxies and their host environ- 
ments. In Section 2 we outline the characteristics of our SDSS sam- 
ple, while in Section 3 we detail our selection of quenched galaxies. 
Readers less interested in the technical aspects of our sample selec- 
tion may wish to begin with Section 4, where we focus on compar- 
ing the properties of quenched galaxies with other populations in 
the SDSS. 

Throughout this paper we adopt a flat cosmology with Q. A = 
0.7, fi M = 0.3 and H = 70 km s -1 Mpc~'. All magnitudes are 
given in the AB system, and we adopt a notation such that equiva- 
lent widths in emission are negative. 



2 DATA 

This work is based on data from the SDSS, supplemented by sev- 
eral value-added galaxy catalogues. We adopt photometric mea- 
surements from the catalogues of Simard et al. (201 1), and galaxy 

' For simplicity we will refer to K+A, E+A and post-starburst galaxies 
collectively as post-starburst galaxies, though we note that this is not an ob- 
servational definition, relying instead on the interpretation of such galaxies' 
spectra as arising from a starburst. 
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stellar masses are estimated from fits to galaxies' ugri spectral en- 
ergy distributions (SEDs) as described by Mendel et al. (2012). 
We adopt measurements of specific star-formation rate (SSFR = 
SFR/M*) from the MPA/JHU data catalogues 2 . In order to provide 
a more robust classification of galaxies' morphology we supple- 
ment the purely quantitative morphological parameters estimated 
by Simard et al. (201 1) with visual classifications from the Galaxy 
Zoo first data release 1 (Lintott et al. 201 1). Below we briefly sum- 
marise the relevant characteristics of these catalogues. 



2.1 Galaxy sample 

The sample used here is selected from the SDSS Data Re- 
lease 7 (DR7; Abazajian et al. 2009) to include all objects with 
extinction-corrected r-band Petrosian magnitudes, m r , satisfying 
14 < m r < 17.77, and classified as galaxies both photometri- 
cally and spectroscopically by the SDSS data reduction pipelines 
(i.e. photoOb j . type = 3 and specObj . specClass = 2). We 
further require that all spectra have a redshift confidence zConf 
> 0.9. Our identification of quenched galaxies requires an estimate 
of the current star-formation rate, and leads us to consider only 
those galaxies with z < 0.2, ensuring that Ho- stays blue-ward of 
significant night-sky emission at A > 8000A. We also adopt a lower 
redshift limit of z = 0.01 to ensure that distance measurements are 
cosmological. 



2.2 Photometric properties 

Simard et al. (2011) performed bulge-disk decompositions for 
galaxies in the SDSS DR7, resulting in a catalogue of structural 
properties for ~ 1 . 12 million objects with extinction-corrected r- 
band Petrosian magnitudes 14 < m r < 18. A full description of 
these decompositions can be found in Simard et al. (201 1); below 
we provide an outline of their methodology to highlight relevant 
details. 

Simard et al. (2011) reprocessed the corrected SDSS im- 
ages (output by the frames pipeline), using sextractor (Bertin & 
Arnouts 1996) to identify objects and generate their associated seg- 
mentation images. Local sky levels were then re-determined on 
an object-by-object basis using a minimum of 20,000 sky pix- 
els, excluding any pixels within 4" of any primary or neighbour- 
ing object mask. Structural decompositions were performed using 
the gim2d software package (Simard et al. 2002) and three differ- 
ent PSF-convolved axisymmetric models: (i) a single-component 
Sersic profile, (ii) a two-component de Vaucouleurs bulge plus ex- 
ponential disk and (iii) a two-component Sersic bulge plus expo- 
nential disk. In all cases fits were performed simultaneously in the 
g and r bands by forcing bulge (half-light size, ellipticity and po- 
sition angle) and disk (scale length, inclination and position angle) 
structural properties to be the same in both bandpasses; remain- 
ing parameters — e.g. total flux, bulge -to-total flux ratio (B/T) and 
model centroid relative to the SDSS object position 4 — were fit in- 
dependently in each bandpass. 



2 accessed at http : //www.mpa-garching .mpg . de/SDSS/DR7/ 

3 accessed at http : //data . galaxyzoo . org. 

4 Simard et al. found that small errors in the image registration between the 
g and r bands resulted in a significant number of galaxies with "positive- 
negative" residuals in their cores. 



In the present work we are primarily interested in galax- 
ies' global properties, and therefore adopt photometric and struc- 
tural parameters from single-component Sersic profile fits. AH 
rest-frame quantities are computed at z = by using k-correct 
v4_2 (Blanton & Roweis 2007) and have been corrected for fore- 
ground (Galactic) extinction using the reddening maps of Schlegel, 
Finkbeiner & Davis (1998). 



2.3 Visual morphologies 

The quantitative morphologies described in Section 2.2 allow for 
analyses of galaxy structural properties, however it is also help- 
ful to relate these quantitative measurements to galaxies' appar- 
ent visual morphology. In general, visual classifications of galaxy 
morphology are carried out by a small number of professionally- 
trained astronomers. While this approach ultimately produces ro- 
bust estimates of galaxy morphology, it is hamstrung by relying 
upon a small pool of identifiers to produce consistent morphologi- 
cal classes across an entire catalogue. 

In order to circumvent this limitation the Galaxy Zoo project 
(Lintott et al. 2008, 201 1) has taken a distributed approach to mor- 
phological classification, leveraging public interest to visually clas- 
sify ~900,000 galaxies from the SDSS. Whereas classical morpho- 
logical typing hinges on a small number of accurate visual types, 
Galaxy Zoo relies on a significant number of less-certain classi- 
fications per galaxy to provide probabilistic estimates of galaxy 
morphology. These data have already been used to explore many 
aspects of the relationship between galaxies and their host environ- 
ments (e.g. Skibba et al. 2009; Bamford et al. 2009; Darg et al. 
2010a,b; Masters et al. 2010). Galaxy Zoo asked users to place 
galaxies into one of four categories: elliptical, spiral, merger and 
"don't know". We use the de-biased elliptical (E) and combined 
spiral (CS) 5 probabilities, which account for the preferential classi- 
fication of galaxies as early type with increasing redshift (see Bam- 
ford et al. 2009 and Lintott et al. 201 1 for details of the de-biasing 
procedure). 



2.4 Stellar masses 

Mendel et al. (2012) provide stellar mass estimates for galaxies in 
the Simard et al. (201 1) catalogue based on fits to their ugri SEDs. 
These fits are conducted using a large grid of synthetic photome- 
try which encompasses a range of plausible age, metallicity, star- 
formation history and dust content. Models are generated using the 
flexible stellar population synthesis (FSPS) code of Conroy, Gunn 
& White (2009), adopting a Chabrier (2003) stellar initial mass 
function (IMF). Estimates of the stellar mass are taken as the me- 
dian of the (marginalized) posterior probability density function, 
with uncertainties quoted as the 16th and 84th percentiles of the 
distribution. Statistical uncertainties on the derived stellar masses 
are of order 0.15 dex, and Mendel et al. show that reasonable varia- 
tions in the 'default' FSPS parameters lead to systematic offsets of 
less than ±0. 1 dex in stellar mass for most galaxies. We adopt here 
masses from the Mendel et al. (2012) catalogue based on the Sersic 
profile fits of Simard et al. (201 1). 



5 Galaxy Zoo asks participants to classify spiral galaxies into sub-classes 
of clockwise, anti-clockwise and edge-on, all three of which are included in 
the combined spiral classification. 
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2.5 Star-formation rates 

Brinchmann et al. (2004, hereafter B04) estimate star-formation 
rates for SDSS galaxies using a combination of detailed spectral 
modelling and empirically-determined calibrations to account for 
contamination from non-stellar ionising sources, e.g. active galac- 
tic nuclei (AGN). The adopted approach depends on the relative 
strength and signal-to-noise (5 IN) of the Hp, [O m]/15007, Ha and 
[N n]/16584 emission lines, which we describe in more detail below. 

Galaxies with well-detected emission lines (S /N > 3) are 
first classified as either star forming, composite or AGN using the 
[Om]^5007/H/3 versus [Nn]i6584/Hff line diagnostic of Baldwin, 
Phillips & Terlevich (1981, hereafter BPT) and the demarcations 
of Kauffmann et al. (2003a) and Kewley et al. (2001); roughly 40 
per cent of galaxies are classifiable in this way. The remaining 60 
per cent are therefore classified based on reduced sets of emission 
lines. Galaxies with S IN > 3 in both Ho- and [Nn]/16584, and 
[Nii]/16584/Hq' > 0.6 (8 percent of the sample) are classified as 
low-S/W AGN. Remaining galaxies with S /N > 2 in Ho- (20 per 
cent) are classified as low-S/W star forming, while galaxies fail- 
ing all of the above criteria are unclassifiable on the basis of their 
emission lines. 

SFRs for star-forming galaxies are determined from detailed 
fitting of their emission lines using the models of Chariot & 
Longhetti (2001). B04 compute the full probability density func- 
tion of a given parameter (e.g. star-formation rate, metallicity, etc.) 
using a Bayesian approach to determine the likelihood of each point 
in a grid of models given the observed data. Although this fitting 
approach incorporates information from multiple emission lines, in 
practice star-formation rates are driven primarily by Ho- luminosity. 

All other galaxies — i.e. those classified as AGN, composite or 
unclassifiable — have SFRs estimated from an empirical relation- 
ship between 4000A break strength (D„4000; Balogh et al. 1999) 
and SSFR based on the star-forming galaxy sample. The advantage 
of this approach is that it allows the determination of star-formation 
rates even for galaxies whose emission-line spectra are heavily con- 
taminated by AGN. However, it is important to note that, since 
D n 4000 is primarily correlated with luminosity-weighted age rather 
than SSFR, adopting the empirical relationship between SSFR and 
D n 4000 of star-forming galaxies automatically assumes that their 
recent star-formation histories are representative of all galaxies 
with a given 4000A break strength. This assumption is particularly 
misleading for galaxies with discontinuous star-formation histo- 
ries, as is the case for so-called 'post-starburst' spectral types (e.g. 
Zabludoff et al. 1996; Goto et al. 2003). We will revisit this issue 
in Section 3. 

For the majority of this work we limit ourselves to using SS- 
FRs determined within the SDSS fibre aperture in order to be con- 
sistent with our own measurements of galaxies' stellar populations 
(Section 3); however, B04 also provide estimates of galaxies' to- 
tal SFR accounting for star formation occurring outside of the 3" 
SDSS fibre aperture. In general, there is good agreement between 
B04's aperture-corrected star-formation rates and estimates of the 
total galaxy star-formation rate obtained from l/V+optical SED fit- 
ting (Salim et al. 2007). 



2.6 Sampling and volume corrections 

In order to discuss absolute quantities we have to correct for the 
visibility of galaxies in our sample, which can be encapsulated by 
two factors: the effect of our photometric selection (Section 2. 1) on 



the sample volume and the impact of our S /N requirements on the 
spectral sampling rate. 

First, galaxies will be observable within some volume V max 
which depends on their SED shape and overall luminosity. Along 
with their stellar mass measurements, Mendel et al. (2012) pro- 
vide estimates of the minimum and maximum observable redshift 
(Zmin.sED and Zmax.sED) f° r each galaxy on the basis of their best-fit 
SED 6 . We use these estimates along with our adopted redshift lim- 
its (0.01 < z < 0.2) to compute the final volume limit, which is 
constrained by 

Zram = max(z mm-S ED,0.01), (1) 

Zmax = min(z max-S ED,0.2) (2) 

and the survey area. In what follows we correct for these volume 
effects by weighting galaxies according to their observable volume 
as V m l x (e.g. Schmidt 1968). 

Second, in order to characterise galaxies' stellar populations 
we require that their spectra have a S /N of at least 10 in the wave- 
length range of interest, 3800-5400A (see Section 3.1.1). For each 
galaxy we estimate S /N based on unmasked pixels used in our 
spectral decomposition and the SDSS error spectra. The detectabil- 
ity of a given galaxy in the SDSS fibre depends on its apparent 
magnitude, concentration 7 and redshift. In order to correct for the 
spectral sampling rate we compute the fraction of galaxies within 
a window of Az = 0.01, AC = 0.2 and Am g = 0.2 that meet our 
S/N requirements, fs/N, and weight galaxies as wssr = /s/Iv- The 
majority of galaxies (94 per cent) have w SSR ^ 2. In the following, 
all data are weighted by wssRKnax to account for sampling effects 
unless otherwise stated. 



3 SELECTING QUENCHED GALAXIES 

Our goal in this work is to study the most recent additions to the 
passive galaxy population; in simple terms, this goal can be met 
by identifying galaxies that are young but host little or no ongoing 
star formation. These criteria can be thought of as analogous to the 
selection of post-starburst galaxies, which is typically based on the 
presence of strong Balmer-line absorption and relative lack of op- 
tical emission lines (e.g. Couch & Sharpies 1987; Zabludoff et al. 
1996; Goto et al. 2003; Blake et al. 2004; Wong et al. 2012), gener- 
alised to identify a more inclusive sample of quenched systems. We 
will discuss a more direct comparison with post-starburst galaxies 
in Section 5.2. 

We first identify galaxies that host little or no ongoing star for- 
mation. As discussed in Section 2.5, the B04 SFRs are only com- 
puted directly from the Ha emission-line flux for galaxies classified 
as star forming on the BPT diagram, or with Ho- S /N > 2 once the 
majority of AGN galaxies are excluded; remaining galaxies have 
their SFR estimated from D„4000. In terms of classifying quenched 
galaxies this is problematic: while Ha emission is sensitive to the 
lifetime of Hn regions, ~10 7 yrs, D„4000 traces the luminosity- 
weighted age of the population, which is sensitive to the integrated 
mass of young stars on timescales over which the stellar mass-to- 
light ratio evolves significantly, i.e. hundreds of Myrs. 

Despite the above complication, pursuing a selection based 

6 Note that the minimum observable redshift arises because the Simard 
et al. (201 1) photometric catalogue excludes galaxies with m r < 14. 

7 C = R9Q/R50, where R90 and R50 are defined as the radii enclosing 90 
and 50 per cent of the r-band Petrosian flux 
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on SSFR has two benefits. First, there is a tight relationship be- 
tween SFR and mass for star-forming galaxies out to at least z = 1 
(e.g. Noeske et al. 2007), which suggests that adopting a division 
between star-forming and passive galaxies based on their instan- 
taneous SFR (or in this case, SSFR) offers a robust selection at 
any redshift. Second, while the canonical selection of post-starburst 
galaxies is based on the absence of particular emission lines, typi- 
cally [O n] or Ha, Yan et al. (2006) have shown that the genesis of 
emission lines, particularly [On], in many post-starburst galaxies 
is not residual star formation but low-ionisation nuclear emission- 
line region (LINER) activity. Selecting quenched galaxies based on 
the absence of [O n] emission may therefore exclude a significant 
fraction of 'true' quenched systems. While a more general selec- 
tion of quenched galaxies can be obtained based on Ha equivalent 
width (Yan et al. 2006), it is important to note that samples selected 
on the absence of Ha still exclude strong-line AGN and almost ex- 
clusively sample the population of LINER-like galaxies (Cid Fer- 
nandes et al. 2010, 2011). In addition, it has been argued that a 
significant fraction of observed LINER-like activity may be due to 
photoionisation by hot evolved stars (e.g. white dwarfs and post- 
asymptotic giant branch [AGB] stars), which could account for the 
extended LINER-like emission observed in many early-type galax- 
ies (e.g. Stasinska et al. 2008; Sarzi et al. 2010; Yan & Blanton 
2012). In this case we might actually expect an excess of LINER 
systems within ~ 1 Gyr of having their star formation quenched, as 
stars in the range 3-6 M Q contribute significantly to the post- AGB 
population (e.g. Taniguchi, Shioya & Murayama 2000; Stasinska 
et al. 2008; Cid Fernandes et al. 201 1). 

We therefore adopt a definition of passive galaxies based on 
both the B04 SSFR and observed EW(Ha) in an effort to be as 
inclusive as possible of potential quenched galaxies. We begin by 
adopting a relationship between fibre SSFR and fibre stellar mass 
for star-forming galaxies in the B04 sample, which is well de- 
scribed by SSFR oc M/ with = -0.16. We identify 329,341 
galaxies with SSFRs a factor of 10 below this fiducial sequence 
(~3cr offset from the SSFR-mass relation), corresponding to a cut 
in D n 4000 of ~1.6. We then supplement this passive sample by in- 
cluding an additional 14,874 galaxies with EW(Ha) > -3A not 
identified as star forming or low-5/A' star forming by B04 8 , cor- 
responding to the EW(Ha) limit described by Cid Fernandes et al. 
(2011) for selecting 'retired' galaxies. 



3.1 Quantifying "young" 

Having outlined the basis for our selection of passive galaxies 
above, we now turn to our classification of galaxies as "young". 
As discussed previously, our method for classifying galaxy stellar 
populations can be thought of as an extension of criteria used to 
select post-starburst galaxies, which typically use cuts on Balmer- 
line equivalent width (EW), e.g. EW(H<5) > 5A, to identify those 
galaxies hosting significant A-star populations. However, adopting 
such strict cuts to identify truncated star-formation histories primar- 
ily selects a subset of quenched galaxies with the highest fraction 



Salim et al. (2007) have shown that there is good agreement between their 
U V-derived star-formation rates and those of B04 for both star-forming and 
low-S/N star-forming galaxies with D„4000 > 1.7, while for D„4000 < 
1.7 B04 may overestimate the SFRs of \ow-S/N star-forming galaxies. As 
our adopted division in SSFR can be effectively translated into a cut in 
D n 4000 > 1.6, any remaining low-5/W galaxies are likely to be genuinely 
star forming. 
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Figure 1. Template spectra and their associated continuum estimates. 
Dashed (red) lines show the broadband continuum estimated as described 
in Section 3.1. Spectra are offset in flux for clarity. 

of young stars and most rapidly declining star-formation rates (e.g. 
Snyder et al. 201 1), and it is unclear how representative these are 
of the majority of galaxies in transition between star forming and 
passive. More recently, various spectral decomposition techniques 
have been used to identify galaxies hosting even relatively minor 
populations of young stars, resulting in a more inclusive sampling 
of the transition population (e.g. Ellingson et al. 2001; Quintero 
et al. 2004; Yan et al. 2006; Nolan, Raychaudhury & Kaban 2007; 
Wild et al. 2009; Yan et al. 2009). It is this latter approach that we 
adopt here. 

We follow Quintero et al. (2004) and Yan et al. (2006) in mod- 
elling galaxy spectra as the simple linear combination of two stellar 
populations, one young and one old. Our underlying old template 
is a 7 Gyr old solar metallicity single stellar population (SSP) com- 
puted using the galaxev stellar population synthesis (SPS) code 
(Bruzual & Chariot 2003, hereafter BC03). The template was cho- 
sen to match the median luminosity-weighted age of galaxies in the 
SDSS (Gallazzi et al. 2005), though our results are relatively insen- 
sitive to the particular choice of SSP template age. A young pop- 
ulation is taken 300 Myr after a starburst lasting 100 Myr, which 
matches the simulated duration of merger-induced starbursts (e.g. 
Snyder et al. 201 1). Our young template is selected to identify the 
strong Balmer absorption associated with A stars, and is therefore 
most sensitive to stellar populations on the order of -0.5-1.5 Gyr 
old. 

Prior to fitting, model spectra are broadened to match the 
SDSS resolution, and all spectra (templates and observed) are lin- 
early rebinned to a dispersion of 1A pixeL 1 . To avoid potential 
mismatches in the continuum shape between observed and tem- 
plate spectra we remove the broadband continuum using a mov- 
ing 200A wide sigma-clipped window. In order to visually match 
the steep spectral shape of our burst template blue-ward of the 
4000A break we adopt a slightly boosted 65th percentile flux (as 
opposed to a true median) for determining the continuum level. The 
resulting continuum is then smoothed using a 100 A boxcar filter; 
examples of the continuum fits to our template spectra are show in 
Figure 1 . Emission lines are masked using an adaptive window that 
scales with the velocity dispersion of the galaxy 9 . Finally, we re- 

9 The standard SDSS spectrold reduction only produces velocity dis- 
persion estimates for galaxies with an early-type spectral classification 
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move any pixels in the observed spectrum with non-zero values in 
the SDSS bad pixel mask: spectra with fewer than 800 unmasked 
pixels (i.e. fewer than half in our wavelength window) are excluded 
from the fitting procedure altogether. 

Decompositions are performed over the wavelength range 
3800-5400A using a non-negative least squares (NNLS) rou- 
tine, and are therefore constrained to physical solutions (i.e. non- 
negative spectral components). The decomposition routine returns 
the weight (by mass) associated with each of the young and 
old templates, w young and w oli respectively; an estimate of the 
mass- weighted young stellar fraction, / yolmg , is therefore given by 
Wyoung /(Wy 0un g + vv id)- Approximately 2 per cent of galaxy spec- 
tra either contain too few unmasked pixels or the NNLS algorithm 
fails to converge on a solution, and these objects are necessarily ex- 
cluded from subsequent discussions. We find that the vast majority 
of galaxies are well described by our simple two-component spec- 
tral model, with typical reduced-^ 2 values of ~ 1.1, and we find no 
evidence for significant variation in the quality of decompositions 
(in terms of x 1 ) with either stellar mass or derived / yoU ng- 

Star-formation histories are known to vary systematically as 
a function of galaxies' stellar mass, such that massive galaxies are 
older and form stars more rapidly than their low-mass counterparts. 
In the context of identifying quenched galaxies, this leads us to 
adopt a relative definition of what makes a galaxy "young" so that 
we can account for these systematic variations in our selection. In 
Figure 2 we show the distribution of / yollng in bins of stellar mass for 
galaxies with 9.5 < log(M/M G ) < 1 1.5. We see that both the mean 
/young and the overall fraction of young galaxies decrease with in- 
creasing stellar mass, as we might expect from the (known) corre- 
lation between stellar mass and age. In each bin we model the / yt) ung 
distribution as the sum of two Gaussians (dashed and dotted lines 
in Figure 2). In order to limit the influence of modelling uncertain- 
ties on our classification we exclude galaxies with extreme values 
of /young (shown as light points in Figure 2). The vertical dashed 
line show the point where the "younger" (i.e., higher / yolmg ) Gaus- 
sian dominates, and above this point we classify galaxies as young. 
In general, our selection of young galaxies is limited to those with 
/young greater than 1 to 2 per cent by mass (i.e. / yollng ~ 0.01-0.02). 

3.1.1 Reliability of / young 

Before proceeding it is useful to assess the overall reliability of 
our / y oung measurements as a function of both spectral S /N and 
/youns- To this end, we have performed a set of Monte Carlo simula- 
tions where we randomly generate composite spectra with a range 
of /young and S/N, and then process these synthetic data using the 
procedures outlined above. We show the results of these simula- 
tions in Figure 3, where the size of individual cells indicates either 
the systematic (top panel) or statistical (bottom panel) uncertainty 
in the measured / yol , ng as a function of simulated (true) / yolmg and 
5 IN. As one might expect there is a strong dependence of the un- 
certainty in /oung as a function of S IN, such that parameters are 
significantly unreliable — both systematically and statistically — at 
low S IN. Relative to the observed distribution of / young (top right 
panel of Figure 3), significant systematic offsets are only avoided 



(eClass < S). As such, we adopt velocity dispersion measurements from 
the Princeton specBS pipeline, accessed at http://das.sdss.org/. In 
general the agreement between spectrold and specBS is excellent. We 
refer the reader to http ://www. sdss . org/dr7/algorithms/veldisp . 
html for details. 



for spectra with S/N > 10, which leads us to adopt this as a S /N 
requirement for the remainder of this work. 

Alternatively, we can ask about the sensitivity of our derived 
results to the parameters used in our spectral decompositions. In 
particular, we would like to know if (and how) our choice of burst 
template age (300 Myr) or metallicity (Z = Z Q ) may bias our over- 
all classifications of galaxies as young. In Figure 4 we plot the de- 
rived fraction of young galaxies as a function of stellar mass, where 
young galaxies are classified as described in Section 3.1. The thick 
solid line shows the fraction of young galaxies as a function of stel- 
lar mass derived using our fiducial templates. For comparison, the 
thin dotted, dashed, dot-dashed and solid lines show the fraction of 
young galaxies determined from fits in which either the burst tem- 
plate age or overall metallicities are varied as indicated in the Fig- 
ure. We find that, while changing the template ages or metallicities 
leads to a factor of ~2 variation in the absolute value of f yomg — that 
is, the inferred mass fraction of young stars — the overall/racfion of 
galaxies classified as young is relatively insensitive to variations in 
the underlying templates due to our relative (as opposed to abso- 
lute) qualification of "young". 

Finally, it is important to note that the fixed 3" aperture of 
the SDSS fibres leads to a variation in physical scale of ~0.6 to 10 
kpc over the redshift range of our sample. We must therefore be 
concerned that any radial variation in galaxies' age or metallicity 
could lead to a systematic bias in our measurements of / young . We 
search for any such dependence in by examining the relationship 
between / y0 ung and z/z milx , where z miix is the maximum redshift at 
which a galaxy could be included in our sample, discussed previ- 
ously in Section 2.6. Our data show evidence for an increase in 
/young with increasing redshift, although any such variation in pas- 
sive galaxies is relatively small (<1 per cent across the volume of 
the sample) and does not appear to significantly affect our spectral 
classifications. The observed increase in / yoU ng is consistent with 
an increasing covering fraction at higher redshift coupled with the 
negative colour gradients observed for many early-type galaxies 
(e.g. Peletier et al. 1990; Goudfrooij et al. 1994; Roche, Bernardi 
& Hyde 2010). While it does not appear that aperture effects sig- 
nificantly influence our results, we nevertheless reiterate that our 
stellar population measurements are restricted to the central few 
kiloparsecs in most galaxies, and hence do not necessarily reflect 
their global properties. 

3.2 Quenched galaxy sample 

Our final sample of quenched galaxies comprises all passive galax- 
ies, as determined by their SSFR, that show evidence for a young 
stellar population, defined as above. In addition, we limit our sam- 
ple to galaxies with log (M/M Q ) > 9.5. All told, we identify 12,105 
such systems, 2.5 per cent of our parent SDSS sample, with a space 
density of -2.2 x 10~ 4 Mpc~ 3 . 

The primary factor governing our selection of quenched galax- 
ies is the rate at which star formation is truncated. In order to illus- 
trate this point, we construct a suite of synthetic star-formation his- 
tories described by 8 Gyrs of constant star formation followed by an 
exponential decline; by varying the characteristic decay rate of the 
decline, t, we can simulate a wide variety of quenching timescales. 
Such models have been shown to provide a relatively good match 
to observed properties of the star-forming galaxy population (e.g. 
Brinchmann et al. 2004). In Figure 5 we plot the timescale over 
which galaxies would be observable given our quenched galaxy 
selection criteria — in this case, the timescale for which a given 
model has / you ng i 0.01 and log SSFR < -11 — as a function of 
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Figure 2. Distribution of /young m 0.1 dex wide bins of stellar mass for galaxies in our spectroscopic sample, where the mean mass is given in the upper 
left of each panel. Open circles show the observed, VJJJjj corrected distribution of /young, while the solid, dashed and dotted curves show the best fitting 
two-component Gaussian. The vertical dashed line indicates where the "younger" (i.e. higher /young) Gaussian dominates, above which we classify galaxies 
as young; see Section 3.1 for details. Light circles indicate points in the tails of the distribution that are excluded from fits. Uncertainties on the observed data 
are computed based on galaxies' individual VJ^ measurements added in quadrature. 



t. This highlights the clear dependence of our sample selection on 
quenching timescale: galaxies whose star formation stops abruptly 
are readily identified as quenched, while galaxies whose star for- 
mation declines over >Gyr timescales are likely not included in 
our quenched galaxy sample. Alternatively, we can cast these sim- 
ulations in terms of their detectability by scaling their observable 
lifetimes according to the total time that quenched galaxies spend 
in transition between the star-forming and passive populations. 
Here we adopt a relatively simple definition for galaxies' transi- 
tion lifetime based on the time that simulated spectra spend with 
1.4 < D„4000 < 1.7, i.e. with intermediate spectral properties. The 
results of this scaling are indicated on the right-hand axis of Figure 



5 by open squares. For quenching timescales of r < 0.5 Gyr we find 
that galaxies are classifiable as quenched for ~ 10 to 60 per cent of 
their transition lifetimes, depending on t. 

The extent to which our sample is representative of the "true" 
transition population therefore depends sensitively on the typical 
timescales over which star formation declines. Kauffmann et al. 
(2004) argue that the lack of a clear environmental dependence 
in the joint distribution of H<5 A and D„4000 points towards slower 
timescales for galaxies' transition from star forming to passive, >1 
Gyr, in agreement with the median timescale of ~ 1.5 Gyrs derived 
by Kaviraj et al. (2007) based on the evolution of UV and optical 
colours. However, a slow evolution for most galaxies seems to be 
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Figure 3. Results of the Monte Carlo simulations described in Section 3.1.1 
for the recovery of /young- In the top left panel we show the systematic 
offset in / yoU ns as a function of input / yoU na anQ S where cell size scales 
with the magnitude of the offset (as indicated) up to a maximum factor of 
2 uncertainty. Open (red) and filled (black) squares indicate negative and 
positive offsets, respectively. For comparison, the top right panel shows the 
observed distribution of / yoU ng ■ In the bottom panel we show the statistical 
uncertainty as a function of f young and S/N, where cells are scaled as above. 
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Figure 4. Fraction of young galaxies as a function of stellar mass for differ- 
ent spectral templates. The thick solid line shows our fiducial decomposi- 
tion, where the burst is 300 Myr old and both components have solar metal- 
licity. Dotted, dot-dashed, dashed and solid lines show the young fractions 
derived for variation in the burst template age of ±200 Myr or the metallic- 
ity of both templates by a factor of 2.5. 




Quenching timescale, x [Gyr] 

Figure 5. Observability of simulated quenched galaxies as a function of 
the quenching rate, characterised by the e-folding time, r. The left axis 
and open (black) circles show the detectability of quenched galaxies given 
our selection criteria (Section 3.2) for values of the decay rate. The right 
axis and open (blue) squares indicate the fraction of their total transition 
lifetime — defined here as the time spent with 1.4 < D n 4000 < 1.7 — for 
which simulated galaxies are detectable as quenched. 



at odds with the persistent bimodal structure of both the SSFR and 
colour distributions: if most galaxies transition from star-forming to 
passive over time scales >1-1.5 Gyrs, then such galaxies should fill 
in the apparent valley at intermediate SSFR and colour (e.g. Balogh 
et al. 2004; Wetzel, Tinker & Conroy 2012). Of course, the most 
likely scenario is that there is no single t that describes galaxies' 
transition from star-forming to passive, and therefore it is more in- 
structive to consider the distribution of r, as in Kaviraj et al. (2007) 
and Martin et al. (2007). Martin et al. (2007) use a combination 
of l/K-optical colour and absorption indices (US a and D„4000) to 
show that the majority of transition galaxies at low redshift (~60- 
80 per cent) are best described by models with t < 0.5-1 Gyr once 
selection effects are accounted for, and that a relative minority of 
galaxies are found with r > 1 Gyr. Based on Figure 5, we would 
therefore expect to identify a majority of quenched galaxies with 
our adopted selection criteria (>50 per cent), with the caveat that 
our results are biased towards those galaxies with the most rapid 
truncation timescales, e.g. r < 0.2-0.3 Gyrs. 

We can obtain a more quantitative picture of our quenched 
galaxy sample by considering their spectroscopic properties rela- 
tive to the general star forming and passive populations. As dis- 
cussed previously, H<5 absorption is primarily sensitive to stars with 
intermediate spectral types (A through early F), and therefore traces 
a particular mass range of stars as they evolve over ~0.5-1.5 Gyr 
timescales. By comparison, the spectral break at 4000 A occurs as 
a result of metal-line blanketing in the atmospheres of cool stars 
blue-ward of 4000A; in hot stars these metal species are multi- 
ply ionised, and therefore the strength of the 4000A break traces 
the overall ageing of a stellar population as progressively cooler 
stars dominate the stellar light. Considering galaxies in the pa- 
rameter space defined by their 4000A-break strength and H<5 ab- 
sorption therefore provides a useful diagnostic for understanding 
the evolution of young stellar populations, and has been used by 
many authors in the past for just this purpose (e.g. Kauffmann et al. 
2003b, 2004; Wild et al. 2007; von der Linden et al. 2010). In Fig- 
ure 6 we plot the distribution of K6 A (Worthey & Ottaviani 1997) 
against a D„4000 for the galaxies in our sample. Open contours 
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Figure 6. The joint distribution of H<5 A and D„4000 for all galaxies in 
our SDSS sample (black contours) and for the spectroscopically-selected 
quenched galaxies (shaded contours). Lines show the predicted evolution 
of quenched galaxies following 8 Gyrs of constant star formation, where 
different colours (orange to blue) indicate different exponential decay rates 
of t = 2, 1, 0.5, 0.2 and 0.05 Gyrs. Open circles show indicate 0.1, 0.5, 1.0, 
1 .5 and 2.0 Gyrs after the onset of quenching. 



show the overall distribution of the SDSS sample, while the shaded 
contours show the parameter space occupied by quenched galaxies. 
The distribution of D„4000 and H(5 A is roughly bimodal, where star- 
forming galaxies occupy a locus with high H<5 A and low D n 4000 
(i.e. the upper-left comer of Figure 6), while old passive galaxies 
have low Hc5 A and high D n 4000. 

The distribution of quenched galaxies in Figure 6 can be un- 
derstood via comparison with the expected behaviour of the trun- 
cated star-formation models discussed previously, which are shown 
in Figure 6 for a variety of e-folding timescales as solid curves. In 
the case of a near-instantaneous truncation of star-formation (or- 
ange line in Figure 6), H(5 A increases as O and B stars terminate 
their evolution and the observed spectrum becomes dominated by A 
stars. The subsequent steady decline of H<5 A and increase of D n 4000 
follows the evolution of the stellar population as it ages. Galaxies 
with the highest values of / youn6 correspond to the strongest Balmer- 
line absorbers (i.e. the upper left of Figure 6), and therefore repre- 
sent the subset of quenched systems caught early in their evolution 
with rapid truncation times (i.e. smallest t). A fraction quenched 
galaxies show extreme values of H<5 A that cannot be accounted for 
in the simple truncated star-formation histories considered here, 
and are more readily described by recent bursts of star-formation; 
these galaxies can be though of as "classical" post-starbursts, and 
make up a relative minority of the sample considered here. 

A more direct comparison between the properties of quenched 
galaxies and the general star-forming and passive galaxy popu- 
lations is shown in Figure 7, where we plot the distributions of 
Hc5 A and D„4000 separately. In both panels, histograms are indi- 
vidually normalised to unity so that the y-axis represents the frac- 
tion of galaxies in a given population. In terms of their D n 4000 
distribution, quenched galaxies form a well defined peak, owing 
to the tight relationship between D„4000 and luminosity-weighted 
age. We also show the distribution of galaxies selected as passive 
based purely on the SSFRs of B04 and those classified as passive 
based on their EW(Ha) as shaded and open (thin) histograms, re- 
spectively. There is a clear division between the two selections at 
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Figure 7. Distribution of D„4000 (top panel) and H<5 A (bottom panel) for 
star-forming, passive and quenched galaxies (dashed, dotted and solid lines, 
respectively). Quenched galaxies are intermediate systems in terms of both 
their D„4000 and H<S A values. The shaded histogram shows the distribution 
of D n 4000 and Hi5 A for quenched galaxies that satisfy our SSFR cut, while 
the thin, open histogram shows the distribution for galaxies with a high 
SSFR from B04, but low EW(Hff) (» -3 A; see Section 3). 



D n 4000 ~ 1.6, which is a direct result of our adopted SSFR limit 
of ~10~ 2 Gyi 1 and the relationship between D„4000 and SSFR 
adopted by B04 (their figure 11). 

The distribution of H6 A in quenched galaxies is generally 
broader than that of D n 4000, in part due to the large uncertainty of 
the H<5 A index relative to D„4000 (cf. figures 5 and 6 of Wild et al. 
2007). Nevertheless, quenched systems occupy a relatively distinct 
peak of US a intermediate between passive and star-forming galax- 
ies. 



4 THE PHYSICAL PROPERTIES OF QUENCHED 
GALAXIES 

The criteria outlined in Section 3 aim to select those galaxies that 
have recently joined the passive population. In this section, we ex- 
amine the colour, morphology and emission-line (AGN) properties 
of quenched galaxies in order to place constraints on their progen- 
itors and identify plausible mechanisms associated with their tran- 
sition. 
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Figure 8. The colour-stellar mass relation (top panel) and colour dis- 
tribution (bottom panel) for galaxies in our spectroscopic sample with 
log (M/Mq) > 9.5. In the top panel, contours show the colour-stellar 
mass relation for our full sample, while shading shows the distribution of 
quenched galaxies. The bottom panel shows the rest-frame g — r colour 
distributions for all ( thin line) and quenched (thick line) galaxies, taking ac- 
count for the slope in the CMR. In both panels the dashed lines bracket the 
region identified as the 'green valley' (see Section 4. 1 for details). 



4.1 Global colours 

As a starting point, we consider the relationship between galaxy 
colour and stellar mass. As discussed in Section 1, the colour- 
mass relation (CMR) is bimodal, with two separate but parallel 
sequences representing quiescent and star-forming galaxies. In the 
upper panel of Figure 8 we show the CMR for galaxies in our full 
sample (black contours); the region of the colour-mass plane oc- 
cupied by quenched galaxies is indicated by shaded contours. A 
significant fraction of quenched galaxies lie intermediate between 
the blue cloud and red sequence — dubbed the 'green valley' — with 
a tendency to be found near the red sequence. This can be seen 
more clearly in the bottom panel of Figure 8, where we plot the 
g—r colour distribution both for galaxies in the full sample (dotted 
histogram) and only quenched galaxies (solid histogram). 

We can gain a better understanding of how galaxies are dis- 
tributed in the colour-mass plane using a simple form for the mass- 
dependent division between red, green and blue galaxies of 



g - r = A + 0.06[log (M/M a ) - 10], 



(3) 



where 0.06 sets the slope of the CMR. We follow Mendez et al. 
(2011) and Lackner & Gunn (2012) in adopting a width for the 
green valley of 0. 1 mags in g — r, centred on the minimum of the 
CMR at log(M/M ) = 10 (A = 0.60; shown by the green lines 
in Figure 8); red and blue galaxies are defined as those above and 
below the green valley, respectively. 

In terms of quenched galaxies 57, 36 and 7 per cent are clas- 
sified as red, green and blue, respectively, compared to 52, 19 and 
29 per cent for the full SDSS sample. While it is relatively easy 
to understand the dearth of blue galaxies in our quenched sample, 
given the prevalence of star formation in the blue cloud population, 
the relatively high abundance of red galaxies may be somewhat 
surprising. To a large extent this is a result of galaxies' non-linear 
transition from blue to red following the truncation of star forma- 
tion, which leads to galaxies evolving rapidly through intermediate 
colours and piling up as they approach the red sequence. It may 
also be that our selection criteria are biased against including the 
youngest (bluest) systems. For example, if the majority of quench- 
ing occurs as a result of AGN activity, then the youngest quenched 
galaxies (with D n 4000 < 1.6) may still have significant Ho- emis- 
sion (e.g. Schawinski et al. 2007), prohibiting their inclusion in our 
sample. However, it is also important to stress that the majority of 
green galaxies (~80 per cent) are classified as star-forming accord- 
ing to the criteria outlined in Section 3, suggesting a clear distinc- 
tion between colour- and spectroscopically-selected samples (see 
also Woo et al. 2012). 

The green valley is generally assumed to host galaxies for a 
short period of time as they transition from star forming to passive 
(e.g. Martin et al. 2007; Mendez et al. 2011). Indeed, as discussed 
earlier, galaxies' colour must evolve rapidly in order to preserve 
the relative under-abundance of green galaxies relative to blue or 
red (e.g. Balogh et al. 2004). We can investigate this scenario ex- 
plicitly with our data by comparing the relative space densities of 
quenched and green galaxies; although our sample is biased to- 
wards the later stages of quenching (i.e. towards red colours), the 
relative abundance of green and quenched galaxies should be com- 
parable if there is no significant evolution in the quenching rate 
on ~Gyr timescales. Quenched galaxies are a factor of ~6 under- 
abundant compared with a colour-selected sample, with a space 
density of 2.2 x 10~ 4 Mpc~ 3 versus 1.5 x 10~ 3 Mpc~ 3 for green 
galaxies. However, as discussed in Section 3.2, our selection crite- 
ria are predominantly sensitive to galaxies experiencing relatively 
rapid quenching. Figure 5 suggests that we may miss over half of 
quenched galaxies with t < 0.5 due to their short observational 
lifetimes, along with a further 50 per cent of galaxies whose star 
formation is quenched over relatively longer timescales. Adopting 
this factor of ~4 suggests that up to 60 per cent of present-day 
green valley galaxies will evolve onto the red sequence over ~Gyr 
timescales, while the remaining ~40 per cent may represent dusty 
star-forming galaxies rather than true transition objects. 



4.2 Visual and quantitative morphology 

Different physical scenarios make different predictions for how 
star formation is quenched; however, based on purely phenomeno- 
logical arguments, the transition of galaxies from star forming to 
passive must, at some level, carry with it an associated transition 
from disk dominated to spheroid dominated in order to jointly 
preserve the morphology-density and star-formation rate-density 
relations. Examining the morphological properties of quenched 
galaxies therefore offers an important constraint on the relative 
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Figure 9. Distribution of spiral (top panel) and elliptical (bottom panel) 
morphological probabilities for different galaxy samples. Star-forming and 
passive galaxies (dashed and dotted histograms, respectively) are defined 
using the division in SSFR discussed in Section 3, while quenched galaxies 
(solid line) shows passive galaxies weighted by P young . Thick (coloured) 
lines show the distribution of morphological probabilities for all galaxies, 
while the grey histograms show the distributions for galaxies with (b/a) < 
0.5 (i < 60°). 



timescales for morphological transformation and star-formation 
truncation. 

In Figure 9 we plot the distributions of Galaxy Zoo ellipti- 
cal and spiral probabilities for star-forming, passive and quenched 
galaxies. This figure shows that quenched galaxies are intermedi- 
ate between star-forming and passive galaxies in terms of their vi- 
sual appearance, with relatively few confidently classified as spiral 
or elliptical. In addition, there is an excess of quenched galaxies 
with low elliptical probability (or, alternatively, high -P sp j ra i) rela- 
tive to other passive galaxies. This excess is consistent with an 
increased fraction of quenched galaxies hosting both bulges and 
disks, supported by the fact that the disparity between passive and 
quenched galaxies is reduced by considering only face-on galaxies 
(b/a < 0.5; thin grey lines in Figure 9). 

We can better characterise any structural differences between 
the star-forming, passive and quenched subpopulations using the 
quantitative morphological decompositions from Simard et al. 
(2011); we show in Figure 10 the distribution of Sersic indices 
for these different subpopulations. The separation between star- 



1 i 1 1 1 1 i 1 1 1 1 i 1 1 1 1 
— — — Star forming 

Passive 
Quenched 




Figure 10. Distribution of Sersic indices for star-forming, passive and 
quenched galaxies (dashed, dotted and solid histograms, respectively). 
Thick (coloured) lines show the distribution of morphological probabili- 
ties for all galaxies, while the grey histograms show the distributions for 
galaxies with (b/a) < 0.5 (i < 60°). 



forming and passive galaxies in terms of their structure is clearly 
represented here, where star-forming galaxies are predominantly 
disky (n < 1.5) and passive galaxies are more centrally concen- 
trated (n > 2). Most strikingly, there is a significant deficit of low 
Sersic index galaxies (n < 1.5) in the quenched population relative 
to star-forming galaxies, suggesting that the hallmark of quenched 
galaxies is their increased central concetration. This confirms our 
expectation based on the analysis of visual classifications in Fig- 
ure 9, i.e. that quenched galaxies are structurally more similar to 
(already) passive galaxies than star-forming ones, and is consistent 
with a picture where the transition of galaxies from star forming to 
passive (or, alternatively, blue to red) is accompanied by a change in 
their morphologies (see also Mendez et al. 201 1 ; Lackner & Gunn 
2012). As noted in the visual classifications, there is an increase in 



the fraction of quenched galaxies with P S] 



:piial 



> p. 



elliptical 



relative to 



the general passive population (~46 vs. 27 per cent). This increase 
is mirrored by a similar rise in the fraction of galaxies with n < 2.5 
(~42 vs. 24 per cent), and may suggest that many quenched galax- 
ies host disks in addition to significant stellar bulges. 



4.3 Incidence of AGN 

The relationship between star-formation and AGN activity remains 
a hotly debated topic. While from a theoretical perspective AGN 
are thought to play a crucial role in regulating star-formation activ- 
ity, observational evidence for a causal link between AGN activity 
and the suppression or cessation of star formation is tenuous. More 
generally, while it is likely that feedback is related in some way to 
star-formation activity, the exact nature of that relationship remains 
poorly understood. In this picture, quenched galaxies offer an at- 
tractive target for the study of feedback processes based on their 
discontinuous star-formation histories. 

In general, optical classification of galaxies as either AGN 
hosts or star forming is based on the comparison of emission-line 
ratios (e.g. Baldwin, Phillips & Terlevich 1981, hereafter BPT) 
with either empirically-derived or theoretically-determined divi- 
sions between different classes. In the top panel of Figure 1 1 we 
show three such demarcations from Stasihska et al. (2006, hereafter 



© 0000 RAS, MNRAS 000, 000-000 



12 Mendel et al. 



S06), Kauffmann et al. (2003a, hereafter K03) and Kewley et al. 
(2001, hereafter KOI) as thick solid, thin solid and dotted lines, re- 
spectively, while the straight dashed line shows the proposed 'opti- 
mal' division of AGN between LINERs and Seyferts determined by 
Cid Fernandes et al. (2010, see also Kewley et al. 2006). Contours 
show the distribution of galaxies in our full sample with S/N > 3 
in each of H/3, [O m]/15007, Ha and [Nn]i6584, while shading in- 
dicates the location of quenched galaxies satisfying the same S/N 
limits. Perhaps unsurprisingly, the majority of quenched galaxies 
with detectable emission lines are classified as AGN (99.3, 98.3 and 
41.8 per cent for S06, K03 and K01 classifications, respectively), 
and the vast majority of these have LINER-like spectra (~60-80 
per cent depending on classification). As discussed in Section 3, 
the abundance of LINER-like emission in quenched galaxies may 
not be surprising. While our selection of galaxies based on the com- 
bination of SSFR and EW(Ha) is more inclusive than, say, a cut on 
just Ha or [O n] (e.g. Yan et al. 2006, 2009), we are unable to prop- 
erly account for the population of strong-line AGN. As an added 
complication, requiring detection in four emission lines limits the 
extent to which weak-lined systems are even classifiable under the 
BPT scheme as both H/3 and [O m]/15007 are fainter and harder to 
detect than Ha and [Nn]/16584. 

As an alternative to the BPT diagram, Cid Fernandes 
et al. (2010) propose a weak-line classification scheme that uses 
EW(Hq-) in place of [O m]A5001/H/3~ the EWHcm2 or WHAN 
(Cid Fernandes et al. 2011) diagram — allowing the classifica- 
tion of a significant number of galaxies with measurable Ha 
and [Nn]/16584 flux, but poor detections of the weaker H/3 and 
[O m]/15007 lines. Cid Fernandes et al. show that requiring reason- 
able S/N in all four of H/3, [O m]i5007, Ho- and [N n]A6584 results 
in a 33 per decrease in sample size overall, and more than a factor 
of 2 decrease in the number of AGN+LINER systems relative to 
classification based on Ha and [Nn]/16584 alone. For our sample, 
requiring detection of only Ha- and [Nn]/16584 nearly doubles the 
number of classifiable emission-line systems. In the bottom panel 
of Figure 1 1, we plot EW(Ha) vs. [Nii]/16584/Ha for galaxies in 
our full sample, where once again quenched galaxies are shown 
as filled contours. Lines represent the so-called 'optimal transla- 
tion' of the classifications shown in the top panel of Figure 1 1 
into the WHAN parameter space. The observed surfeit of LINER- 
like systems suggested by the BPT classification is confirmed by 
the WHAN classification scheme, where the majority of systems 
are found to have relatively low EW(Ha). However, as cautioned 
above, our ability to properly classify quenched Seyfert galaxies is 
hindered by estimates for the star-formation rates in these galaxies, 
which may be systematically overestimated by adopting D n 4000 
as a proxy (e.g. Salim et al. 2007); our data therefore represents a 
lower limit on the true Seyfert fraction in such systems. 



5 DISCUSSION 

We have used a combination of Galaxy Zoo visual morpholo- 
gies and quantitative structural parameters to study the properties 
of quenched galaxies in the SDSS. Our main result is that these 
quenched systems appear to host significant bulge components de- 
spite their young ages. 

5.1 The relationship between quenching and morphology 

The buildup of the present-day red sequence can be understood 
as a combination of quenching processes and merging (e.g. Faber 
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Figure 11. Line ratio diagnostics for galaxies in our SDSS sample. Top 
panel: BPT diagram for the 143,727 galaxies (~30 per cent) with S/N > 3 
in each of H/3, [O m]/15007, Ha and [Nn]/16584 shown as open contours, 
with quenched galaxies shown as shaded contours. Lines show differ- 
ent divisions between star-forming and AGN galaxies proposed by S06 
(thick solid line), K03 (thin solid line) and K01 (dotted line). The dashed 
line shows the division between Seyfert and LINER galaxies proposed 
by Cid Fernandes et al. (2010). Bottom panel: The equivalent width Ha 
vs. [Nn]/(6584/Ha (WHAN) digram proposed by Cid Fernandes et al. 
(2010) for the 263,530 galaxies (-54 per cent) with S/N > 2 in Ha 
and [Nn]/16584. Lines are the same as in the above panel, but now rep- 
resent the 'optimal transposition' of the K03, S06 and K01 classifications 
into EW(Ha)-[N n]/16584/Ha parameter space as derived by Cid Fernandes 
et al. (2010). 



et al. 2007). In this simplified picture, galaxies' recent merger his- 
tory dictates their morphology — major mergers produce elliptical 
galaxies and minor mergers build spheroids without destroying 
stellar disks — while star formation is governed by a combination 
of feedback and environmental mechanisms. 

In this context, quenched galaxies provide leverage in disen- 
tangling the relative interplay between processes that govern galax- 
ies' supply of cold gas and those that can alter their structure. As 
discussed in Section 4.2, Figures 9 and 10 show that quenched 
galaxies are both visually and structurally similar to other quies- 
cent galaxies, despite their apparent recent arrival into the passive 
population. Post-starburst galaxies, which represent an extreme tail 
of the quenched galaxy spectrum, paint a similar picture: Quintero 
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et al. (2004) show that post-starburst galaxies are more centrally 
concentrated (have higher Sersic indices) than star-forming galax- 
ies, and Balogh et al, (2005) show that while post-starburst galax- 
ies host prominent bulges, their apparent progenitors — Balmer- 
strong emission galaxies, so-called e(a) systems (e.g. Dressier et al. 
1999) — are still largely disk dominated. More generally, Simard 
et al. (2009) showed that, at least in clusters, there is a correlation 
between the fraction of galaxies with [O n] emission (i.e. star for- 
mation) and the fraction of galaxies with late-type morphologies. 

It remains unclear, however, if the observed relationship be- 
tween morphology and star-formation is causal — i.e. that a stellar 
bulge is required for quenching to occur — or if the buildup of a 
stellar bulge is the result of physical processes that also induce qui- 
escence. Given the young ages of galaxies considered here, our re- 
sults suggest that bulge formation must precede the shutdown of 
star formation, i.e. that morphology is a critical factor in determin- 
ing quiescence. Bell et al. (2012) have shown that there is a relation- 
ship between star formation and structure since at least z ~ 2, where 
the buildup of a passive population is mirrored by the buildup of 
galaxies with high Sersic indices. Similarly, Wuyts et al. (2011) 
showed that at any given epoch there is a well defined relationship 
between galaxies' structure and their position on the star forma- 
tion rate-stellar mass plane. However, there is significant scatter 
between quiescence and any number of galaxy properties, suggest- 
ing that the mapping between structure and star formation is, at 
best, imperfect. Indeed, there are notable exceptions to a paradigm 
where galaxies' star-formation activity is directly tied to their struc- 
ture; passive spirals in particular pose a problem to any theory of 
galaxy evolution that limits the formation of passive galaxies to 
violent events, e.g. major mergers, and there is some evidence sug- 
gesting that quiescent disks may be quite common, particularly at 
z > 0.8 (Bundy et al. 2010). However, both Masters et al. (2010) 
and Bundy et al. (2010) show that even passive disk galaxies appear 
to host more prominent bulges than their star forming counterparts, 
consistent with the results presented here. 

5.2 Where do post-starburst galaxies fit in? 

As mentioned previously, our selection of quenched systems builds 
upon techniques primarily used to identify post-starburst galax- 
ies; it is therefore worthwhile to revisit the comparison between 
these two populations in light of the results discussed in Section 4. 
By construction, the stringent criteria adopted for selecting post- 
starburst galaxies — e.g. EW(HS A ) > 4-5A and no [On] or Ha 
emission — identify galaxies with extreme star-formation histories. 
Of the 12,105 quenched galaxies in our sample, only 721 have Hi5a 
equivalent widths > 4A, or ~0. 1 per cent of the total SDSS sample, 
in relatively good agreement with other studies of post-starburst 
galaxies (Goto et al. 2003; Blake et al. 2004; Goto 2005; Wong 
etal. 2012). 

Several studies have highlighted the significant burst-mass 
fractions required to reproduce the absorption-line properties of lo- 
cal post-starburst galaxies. Based on fits to their t/V+optical SEDs, 
Kaviraj et al. (2007) find that 20-50 per cent of post-starburst 
galaxies' total stellar mass may be formed in recent bursts, while 
Norton et al. (2001) and Swinbank et al. (2012) derive more modest 
mass fractions of order 5 to 20 per cent. By comparison, the inferred 
A-star mass fractions from our own fits are substantially lower, of 
order a few per cent, for the the majority of quenched galaxies. In 
addition, we find a minority of galaxies with / yolm e greater than 10 
per cent, and virtually none with / yollI1£ , greater than 20 per cent. 
It could be, then, that post-starburst galaxies exclusively trace very 



efficient star-formation activity — as may be the case in gas-rich ma- 
jor mergers (e.g. Cox et al. 2006; Saintonge et al. 2011) — while 
the broader population of quenched galaxies host comparatively 
low-levels of star-formation prior to becoming quiescent. In fact, 
merger-driven formation is also suggested by the high incidence of 
disturbed morphologies among post-starburst galaxies (Yang et al. 
2004a; Blake et al. 2004; Pracy et al. 2009); the lack of a clear dis- 
tinction between the local environments of post-starburst and 'nor- 
mal' galaxies lends further support to merging over more specific 
environmental process, e.g. ram-pressure stripping (Zabludoff et al. 
1996; Blake et al. 2004; Hogg et al. 2006; Wong et al. 2012). 

Yet, despite any differences between quenched and post- 
starburst galaxies in terms of their formation processes, the end 
products appear remarkably similar. As discussed in Section 5.1, 
one of the hallmarks of quenching appears to be the presence of a 
stellar bulge, and this holds equally true for post-starburst galaxies 
as it does for the more general quenched population (Quintero et al. 
2004; Balogh et al. 2005; Goto 2005; Hogg et al. 2006; Pracy et al. 
2009). However, among quenched galaxies there is a clear corre- 
lation between H6a absorption-line width and central stellar mass 
surface density, pointing to increasingly efficient gas dissipation in 
the most extreme Balmer-line absorbers. It may be, then, that true 
distinction between HS-strong post-starburst galaxies and the re- 
mainder of the quenched galaxy population lies in the efficiency of 
bulge formation, rather than a particular formation mechanism. 

5.3 What quenches star formation? 

Ultimately, we would like to identify the process(es) that drive 
galaxies from star-forming to passive. Given the relatively short 
lifetime of the quenching phase, ram-pressure stripping is an ob- 
vious culprit to account for the rapid decline of star-formation ac- 
tivity in quenched galaxies; however, it has been shown that simple 
gas stripping prescriptions dramatically overestimate the fraction of 
passive satellite galaxies in semi-analytic models (e.g. Weinmann 
et al. 2006b; Kimm et al. 2009). Furthermore, such simple models 
are inconsistent with the more gradual decline in star-formation ac- 
tivity inferred from statistical analyses of galaxy stellar populations 
(e.g. Kauffmann et al. 2004). Balogh et al. (2009) show that adopt- 
ing a slower "strangulation" model for star formation as galaxies 
fall into dark matter haloes reduces some of the tension between 
observed and predicted colour distributions, but results in an over- 
production of green galaxies. Recently, Wetzel, Tinker & Conroy 
(2012) have argued that the invariance of central and satellite SSFR 
distributions points towards a hybrid quenching model, where satel- 
lite galaxies evolve as normal for several Gyrs following accretion 
onto a massive dark matter halo and are ultimately quenched as a 
result of mergers, tidal interactions or harassment. In fact, our data 
are consistent with any number of environmental mechanisms that 
may limit galaxies' gas supply without immediately truncating cen- 
tral star-formation activity; however, differentiating between such 
mechanisms requires a more detailed study of galaxies' environ- 
ments that is beyond the scope of the present work. Nevertheless, 
while the extent to which environment plays an important role in 
the evolution of "average" galaxies remains uncertain (e.g. Peng 
et al. 2010), the scarcity of pure disks in our sample (Figure 10) 
suggests that even if all passive galaxies ultimately form as the re- 
sult of gas stripping processes in high-density environments, ad- 
ditional mechanisms above and beyond simple, wholesale gas re- 
moval are required to account for the correlation between bulge 
formation and quenching (see also Bundy et al. 2010; Bell et al. 
2012). 
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Alternatively, galaxies' star formation may be regulated inter- 
nally by energetic feedback from a central black hole (e.g. Croton 
et al. 2006; Hopkins et al. 2008b,a). As discussed in Sections 3 and 
4.3, our dependence on optical emission-line data to identify star- 
forming galaxies limits the extent to which we can study the as- 
sociation between optical AGN and quenching; however, it is still 
interesting to consider the relationship between quenched galaxies 
and other sub-samples with the understanding that the AGN frac- 
tion in quenched galaxies is a lower limit. 

We find that the fraction of quenched galaxies classified as 
AGN (i.e. having S /N > 3 in H/3, [Om]^5007, Ha and [Nii]A6584 
as well as satisfying the K03 classification line) is comparable to 
that of other "old" galaxies'": 8 per cent of quenched galaxies are 
classified as AGN, compared to 17 and 7 per cent of young and 
old galaxies, respectively. As discussed in Section 4.3, the ma- 
jority of AGN in quenched galaxies appear LINER-like (70-80 
per cent) and have L[omj < 10 7 L Q ; in contrast, AGN in young 
galaxies are more than twice as likely to have high [O m] luminosi- 
ties (>10 7 L ) than either quenched or old systems. These results 
are in good agreement with previous studies showing that AGN 
are more frequently associated with recent or ongoing star forma- 
tion (e.g. Kauffmann et al. 2003a; Schawinski et al. 2007; Best & 
Heckman 2012). We obtain a similar picture if we consider the fre- 
quency of high- and low-excitation radio-loud AGN using the re- 
cent catalogue of Best & Heckman (2012), which suggests that the 
paucity of AGN in quenched galaxies may be genuine (see also 
Shin, Strauss & Tojeiro 2011). It is interesting that we find no ev- 
idence for increase in the AGN activity of quenched galaxies, and 
hints that nuclear feedback may not be directly responsible for their 
transition from star forming to passive, however given the limita- 
tions of our sample we cannot make a more conclusive statement. 

Finally, Martig et al. (2009) show that the addition of a 
dynamically-hot stellar component can suppress fragmentation in 
gas disks, so-called "morphological quenching". In this scenario, 
the addition of a massive bulge results in a deepening of the cen- 
tral potential which, in turn, drives an increase in the epicyclic fre- 
quency of the disk; morphologically-quenched galaxies are there- 
fore quenched despite hosting significant gas disks. Morphological 
quenching offers a mechanism to shutdown star formation in the ab- 
sence of environmental and feedback effects, and can account for 
both the significant population of disks in our sample and the ap- 
parent lack of nuclear activity. Such a mechanism may also account 
for the significant population of local early-type galaxies observed 
to host molecular gas (22 per cent; Young et al. 201 1) but little on- 
going star formation. If morphological quenching does contribute 
significantly to the suppression of star-formation in average galax- 
ies, it remains important to understand why bulges appear to be a 
necessary but insufficient component of the quenching process (e.g. 
Bell et al. 2012; see also discussion in Bundy et al. 2010). 



6 CONCLUSIONS 

We have used the SDSS to identify the rare population of galaxies 
in transition between star forming and passive, and have used 
their photometric properties to study the mechanisms that govern 
star-formation quenching. Our results can be summarised as 
follows: 



Here, the classification of galaxies as young or old is based on the dis- 
tribution of firoungi i- e - Figure 2 



i) Quenched galaxies are able to account for up to 60-70 per 
cent of the present day green-valley population after correcting 
for selection effects, suggesting that the majority of green 
galaxies will transition from star-forming to passive on ~Gyr 
timescales. However, the two samples select fundamentally 
different galaxy populations: while we require that quenched 
galaxies are passive, up to 80 per cent of "green" galaxies host 
significant emission. 

ii) Both visually and quantitatively, quenched galaxies appear 
to host significant bulge components despite their young ages. 
It appears that the transition from star forming to passive is 
accompanied by significant morphological transformation, but 
that this transformation takes place in advance of or coincident 
with the shutdown of star formation. 

iii) A comparison between quenched galaxies and classical post- 
starbursts shows that both populations conform to a scenario 
where bulge growth is an integral part of the quenching process. 
Based on the comparison of their central stellar surface densities 
we suggest that the defining characteristic of post-starburst 
galaxies is the efficiency of gas dissipation in the formation of 
their bulges, rather than a particular formation mechanism. 

iv) Limitations of our sample selection notwithstanding, we 
find no evidence for an increase of the optical AGN fraction in 
quenched galaxies relative to other galaxies in the SDSS. These 
results suggest that, if AGN do play a decisive role in quenching 
star formation, then nuclear activity is relatively short lived 
compared to the observable lifetime of quenched galaxies. 

Overall, our results show that galaxy quenching is closely 
tied to morphology, either through processes directly related to the 
presence of a bulge component (e.g. Martig et al. 2009), or indi- 
rectly through associated correlations with, e.g., black hole mass 
(e.g. Hopkins et al. 2008b). With our current data we are unable 
to identify the catalyst for bulge growth — which may be triggered 
internally via secular processes, or externally as the result of tidal 
interactions — however our results suggest that this process is a crit- 
ical component governing the buildup of passive stellar mass with 
time. Future multi-wavelength observations should help to distin- 
guish between the multiple evolutionary pathways currently avail- 
able to galaxies in theoretical models. 
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